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We apply the GiBUU model to questions relevant for current and future long-baseline neu- 
trino experiments, we address in particular the relevance of charged- current reactions for neutrino- 
disappearance experiments. A correct identification of charged-current quasielastic (CCQE) events 
— which is the signal channel in oscillation experiments — is relevant for the neutrino energy re- 
construction and thus for the oscillation result. We show that about 20% of the quasielastic cross 
section is misidentified in present-day experiments and has to be corrected for by means of event 
generators. Furthermore, we show that a significant part of Itt"'" (> 40%) events is misidentified as 
CCQE events mainly caused by pion absorption in the nucleus. We also discuss the dependence of 
both of these numbers on experimental detection thresholds. We further investigate the influence 
of final-state interactions on the neutrino energy reconstruction. 



I. INTRODUCTION 

Good knowledge of the neutrino energy is required for a 
precise determination of oscillation parameters in Vfj, dis- 
appearance measurements. However, the neutrino beam 
is far from being mono-energetic in present experiments. 
Thus, t'^j disappearance experiments search for a distor- 
tion in the neutrino flux in the detector positioned far 
away from the source. Comparing both, un-oscillated 
and oscillated flux, one gains information about the os- 
cillation probability and with that about mixing angles 
and mass squared differences. 

The neutrino energy is not measurable directly but has 
to be reconstructed from the reaction products.^ Present 
oscillation experiments use the charged-current quasielas- 
tic (CCQE) reaction as signal event and reconstruct the 
energy with quasifree two-body kinematics from the out- 
going muon assuming the target nuclcon is at rest. Two 
immediate questions arise from this procedure: (1) How 
good is the identification of CCQE events? (2) How ex- 
act is the assumption of quasifree two-body kinematics 
for nucleons bound in a nucleus where many in-medium 
modifications arc present? 

CCQE is defined as lyn — > £^p (i.e., on the nucleon). 
In the nucleus, CCQE is masked by final-state interac- 
tions (FSI). Thus, the correct identification of CCQE is 
immediately related to the question of how FSI influence 
the event selection. The main background to CCQE is 
CCl7r+ production. If the pion is absorbed in the nu- 
cleus and/or not seen in the detector, these events can 
be misidentifled as CCQE events. Consequently, a proper 
understanding of both CCQE and CCl7r+ on nuclei is es- 
sential for the reconstruction of the neutrino energy. 

This article addresses the questions outlined above in 
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^ External experiments are performed to help predict the beam 

profile. See, for example, the hadron production experiment at 

CERN (HARP) Q. 



the framework of the GiBUU transport model. After giv- 
ing the necessary model details, we start with a general 
introduction of event classification in typical neutrino de- 
tectors and discuss how possible detection thresholds in- 
fluence the measured spectra. We further discuss both 
CCQE and CCItt"'' cross sections and their entanglement. 
Finally, we investigate how nuclear effects influence the 
reconstruction of the neutrino energy. 

II. NEUTRINO SCATTERING IN THE GIBUU 
TRANSPORT MODEL 

The presence of in-medium modifications and, in par- 
ticular, final-state interactions inside the target nucleus 
requires the use of state-of-the-art theoretical methods 
for the extraction of elementary processes from exper- 
iments with nuclear targets (see also the review by 
Alvarez-Ruso [3]). For this aim we use the GiBUU trans- 
port model 0, where the neutrino first interacts with a 
bound nucleon. It has recently been questioned how good 
this impulse approximation (lA) actually is 0, • Our 
experience with photonuclear processes Q indicates that 
this approximation gives reasonably reliable results for 
momentum transfers > 0.2 GeV. The use of the lA re- 
quires a good description of both the elementary vertex 
and the in-medium-modifications. The final state of this 
initial reaction undergoes complex hadronic final-state 
interactions. In the following, we give a brief overview of 
our model and refer the reader in particular to Ref. Q 
and to our website Q for more details. 

In the energy region of i?^ ^ 0.5 — 2 GeV, the elemen- 
tary neutrino-nucleon cross section contains contribu- 
tions from quasielastic scattering (QE: £N — > £'N'), reso- 
nance excitation (R: £N — )• £'R) and direct (i.e., nonreso- 
nant) single-pion production (BG: £N — i'lrN') treated 
in our description as background. QE scattering is the 
most important process at these energies, followed by 
single-pion production through the excitation and sub- 
sequent decay of the A resonance P33(1232). However, 
we include in addition 12 N* and A resonances with in- 
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variant masses less than 2 GeV. The vector parts of the 
single contributions are obtained from recent analyses of 
electron-scattering cross sections. The axial couplings are 
obtained from PCAC (partial conservation of the axial 
current), and, wherever possible, we use neutrino-nucleon 
scattering data as input. 

The elementary neutrino-nucleon cross section is mod- 
ified in the nuclear medium. Bound nucleons are treated 
within a local Thomas-Fermi approximation which natu- 
rally includes Pauli blocking. The nucleons are bound in 
a mean- field potential depending on density and momen- 
tum which we account for by evaluating the above cross 
sections with full in-medium kinematics. We further con- 
sider the collisional broadening of the final-state parti- 
cles within the low-density approximation Fcoii = p<yv 
obtained in a consistent way from the GiBUU cross sec- 
tions. Details of our model for the elementary vertex and 
the corresponding medium modifications can be found in 
Ref. 0. 

In the next step, the particles propagate through the 
nucleus undergoing final-state interactions (FSI) which 
are simulated with the coupled-channel semiclassical 
GiBUU transport model (more information and code 
download on our website [1]). The GiBUU model is 
based on well-founded theoretical ingredients and has 
been tested in various very different nuclear reactions; in 
particular, against electron- and photon-scattering data 

The space-time evolution of a many-particle system in 
a mean-field potential is described by the BUU equation. 
For particles of species i, it is given by 

{dt + VpH ■ Vr - VrH ■ Vp) /,(r,p,t) = 

^oll[/i, /w, /tt, /a, •■•], (1) 

where the phase-space density fi{r,p, t) depends on time 
t, coordinates r and the four-momentum p. H is the rel- 
ativistic Hamiltonian of a particle of mass M in a scalar 
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potential U given hy H = f[M + t/(r,p)] +p2j . 

The scalar potential U usually depends both on four mo- 
mentum and on the nuclear density. The BUU equa- 
tions are coupled through the collision term IcoW which 
accounts for changes (gain and loss) in the phase-space 
density due to elastic and inelastic collisions between par- 
ticles and also to particles decaying into other hadrons. 
In particular, we include two-body reactions like e.g. 
TTiV ^ ttN, NN NN, RN NN, RN R'N, and 
three-body processes like irNN -> NN and ANN — >• 
NNN. By this coupled-channel treatment we can de- 
scribe side-feeding processes into different channels. This 
complex set of coupled differential-integral equations is 
then solved numerically with the GiBUU code. 

All particles (also resonances) are propagated in mean- 
field potentials according to the BUU equations. Those 
states acquire medium-modified spectral functions (nu- 
cleons and resonances) and are propagated off shell. The 
medium modification of the spectral function is based 



both on collisional broadening and on the mean-field po- 
tentials, both of which depend on particle kinematics as 
well as on nuclear density. 

Altogether, FSI lead to absorption, charge exchange 
and redistribution of energy and momentum, as well 
as to the production of new particles. We have shown 
in earlier works, that their impact on neutrino-induced 
pion production is dramatic [1, Thus, a qualitatively 
and quantitatively correct treatment of these effects is of 
great importance, especially for the energy reconstruc- 
tion as we will demonstrate in the following. 

III. EVENT SELECTION 

Event selection in current neutrino experiments is a 
highly complicated subject. Rather than presenting a 
quantitative discussion for each particular setup we give a 
qualitative picture of how nuclear effects themselves mod- 
ify the measured spectra in chargcd-currcnt (CC) scat- 
tering assuming certain detection methods. Two generic 
detectors with the following properties are used toward 
this aim. We note that the event identifications used here 
are the ones used in the actual experiments. 

Cherenkov detector. In a Cherenkov detector (e.g., 
MiniBooNE and K2K-lkt), CCQE events are iden- 
tified by a single ring from the outgoing lepton. 
Muons can be tagged by their decay electron. If 
pions are produced, they lead to additional rings 
either from the 7 decay of the 7r° or from the decay 
muon of the charged pions. 

For the Cherenkov detector, we identify the two 
relevant processes in the following way: 

CCQE: Ifi- 07r+ Ott" 07r° xp xn , 
CCl7r+: Ifi- 1tt+ Ott" 07r° xp xn , 

where xp and xn indicate, that any number of pro- 
tons or neutrons are allowed. 

The lower momentum thresholds depend on the in- 
dex of refraction n via 

-,1,1 m 

Pthrcs = - ^ P thros = , ^ , (2) 

n - 1 

where m is the particle mass. From this, one eas- 
ily obtains the kinetic energy thresholds. Typical 
values for water (n = 1.33) are Tthrcs ~ 55 McV 
for muons, 75 MeV for charged pions, MeV for 
neutral pions (identified via their 7 decay), and 485 
MeV for protons. Lower thresholds (« 10 MeV for 
muons and charged pions, « 65 MeV for protons) 
are reached with the MiniBooNE detector which is 
filled with mineral oil with n = 1.47 [l3l and, in 
addition, produces scintillation light. 

Tracking detector. In a tracking detector (e.g., Sci- 
BooNE and K2K SciFi), all charged particles leave 
tracks which can be used to identify the particles 
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and determine their properties. Thus, highly ad- 
vanced event selection procedures are applied. To 
keep it simple, we identify 

CCQE: 1^- 07r+ Ott" On° Ip xn , 
CCl7r+: 1tt+ Ott" Ott" xp xn . 

The thresholds depend strongly on the experimen- 
tal setup, e.g., the SciFi detector requires both 
muon (pion) kinetic energy to be above « 500 MeV 
(100 MeV) and the proton kinetic energy above 175 
MeV [lil- 

In the following, we assume perfect particle identifica- 
tion above threshold in both cases and neglect any other 
experimental restrictions. 

IV. TOPOLOGIES 
A. CCQE identification 

The CCQE reaction, f^n £~p, being the dominant 
cross section at low energies, is commonly used to recon- 
struct the neutrino energy. In other words, CCQE is the 
signal event in the present oscillation experiments. 

The experimental challenge is to identify true CCQE 
events in the detector, namely, muons originating from an 
initial QE process. To be more precise, true CCQE corre- 
sponds to the inclusive CCQE cross section including all 
medium effects or, in other words, the CCQE cross sec- 
tion before FSI. The difficulty comes from the fact that 
the true CCQE events are masked by FSI in a detector 
built from nuclei. The FSI lead to misidentified events, 
e.g., an initial A whose decay pion is absorbed or which 
undergoes "pion-less decay" contributes to knock-out nu- 
clcons and can thus be counted as a CCQE event — we 
call this type of background event a "fake CCQE" event. 
We denote every event which looks like a CCQE event as 
"CCQE-like". 

As outlined above, in Cherenkov detectors CCQE-like 
events are all those where no pion is detected, whereas 
in tracking detectors, CCQE-like events are those where 
a single proton track is visible and at the same time 
no pions are detected. The two methods are compared 
in Fig. [TJ The "true CCQE" events are denoted with 
solid lines, the CCQE-like events by dashed lines. The 
Cherenkov detector is able to detect almost all true 
CCQE events (top panel; solid vs. dash-dotted lines ap- 
proximately agree) but sees also a considerable amount 
of "fake CCQE" (or "non-CCQE") events (top panel; 
the dashed line is roughly 20% higher than the solid 
line). They are caused mainly by initial A excitation 
as described in the previous paragraph (absorption of 
decay pion or "pion-less decay"); their contribution to 
the cross section is given by the dotted lines. These 
additional (fake) events have to be removed from the 
measured event rates by means of event generators, if 
one is interested only in the true QE events. It is ob- 
vious that this removal is better the more realistic the 
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FIG. 1: (Color online) Total QE cross section on ^^C (solid 
lines) compared to different methods on how to identify 
CCQE-like events in experiments (dashed lines). The top 
panel shows the method commonly applied in Cherenkov de- 
tectors; the lower panel shows the tracking-detector method 
as described in the text. The contributions to the CCQE-like 
events are also classified [CCQE-like from initial QE (dash- 
dotted) and from initial A (dotted lines)]. Experimental de- 
tection thresholds are not taken into account. 



generator is in handling the in-medium ttA^A dynamics. 
On the contrary, less CCQE-like than true CCQE events 
are detected using the method applied in tracking de- 
tectors, which triggers both on pions and protons (lower 
panel, difference between dashed and solid line). The 
FSI of the initial proton lead to secondary protons or, 
via charge exchange to neutrons which are then not de- 
tected as CCQE-like any more {single proton track). We 
find that at tracking detectors the amount of fake events 
in the CCQE-like sample is less than at Cherenkov de- 
tectors (dashed and dash-dotted lines almost agree with 
each other in the lower panel but not in the top panel). 
We conclude that, even if the additional cut on the pro- 
ton helps to restrict the background, an error of about 
20% remains since the measured CCQE cross section un- 
derestimates the true one by that amount. Note that 
experimental detection thresholds are not yet taken into 
account. Thus, about 20% of the total cross section has 
to be reconstructed by using event generators. In this 
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FIG. 2: (Color online) Ratio of the CCQE-like to the true 
CCQE cross section as a function of the lower proton (pion) 
kinetic energy detection threshold for CC i'^ on ^^C at Ei, — 
1 GeV. The solid lines are obtained using the tracking detec- 
tor identification, whereas the dashed lines are for Cherenkov 
detectors. 

case, these generators have to be very realistic in de- 
scribing the in-medium nucleon-nucleon interactions. 

To investigate further the relationship between the 
CCQE-like and true CCQE cross section, we show their 
ratio as a function of the lower proton and pion-kinetic 
energy detection thresholds in Fig. [5] (see Sec. Illll for the 
thresholds applied in present experiments). As the pro- 
ton is not at all relevant for the CCQE identification in 
Cherenkov detectors, the ratio is independent of the pro- 
ton kinetic energy detection threshold (dashed line in top 
panel). This is very different in tracking detectors which 
rely on the detected proton — here the efficiency is re- 
duced to «10% at a proton kinetic energy threshold of 
0.5 GeV (solid line in top panel). Even at T^^^^^ = 0, 
the efficiency does not exceed 80% because of charge- 
exchange processes that lead to the emission of unde- 
tected neutrons and because of secondary proton knock- 
out that leads to multiple-proton tracks. These effects 
cause the difference between the solid and the dashed 
lines in the top panel of Fig. [TJ Focussing on the lower 
panel of Fig. [51 we find that the CCQE-like cross section 
increases for both detector types as T^y^^.^^ increases. In 



FIG. 3: (Color online) Ratio of the QE-induced CCQE-like 
(i.e., the fraction of the CCQE-like event where the initial 
interaction was QE: i/^n — >■ iJ,~p) to the true CCQE cross sec- 
tion as a function of the lower proton (pion) kinetic energy 
detection threshold for CC on ^^C at £^ = 1 GeV. The 
solid lines are obtained using the tracking detector identifi- 
cation, whereas the dashed lines are for Cherenkov detectors. 



this case, even more events with pious in the final state 
appear as CCQE-like because then these pious are below 
threshold and thus not detected. 

The CCQE-like cross section is split into QE and non- 
QE sources (like A excitation) in Fig. [3] and Fig. H) The 
top panel of Fig.[3]shows again the ability of a Cherenkov- 
like detector to identify over 98% of the initial CCQE 
events (dashed line); the missing strength is mainly lost 
into pion channels, that is, the nucleons rescatter and 
produce pions such that the event is no longer classi- 
fied as CCQE-like. This fraction almost vanishes (the 
dashed line gets even closer to one in the lower panel) 
when the pion kinetic energy threshold increases because 
then the CCQE-induced pions are no longer detected and 
the event counts again as CCQE-like. 

Let us now turn to the non-QE CCQE-like cross sec- 
tion displayed in Fig. |4l Different sources are indicated: 
initial A excitation and initial single-pion background re- 
action (higher resonances are negligible here and thus not 
shown). The top panel shows again the dependence on 
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FIG. 4: (Color online) Ratio of the non-QE induced CCQE- 
like (i.e., the initial interaction was not Uf^n -¥ ^.~p) to 
the true CCQE cross section as a function of the proton 
(pion) detection kinetic energy threshold for CC i^^ on ^^C 
at = 1 GeV. The solid lines are obtained using the 
tracking detector identification, while the dashed lines are for 
Cherenkov detectors. The bullets symbolize the cross sections 
where the initial reaction is A excitation, whereas the open 
symbols stand for the initial single-pion background. The 
contribution from higher resonances is negligible. 



the proton-kinetic energy threshold, which is not rele- 
vant in the Cherenkov case, where the non-QE CCQE- 
like contribution adds up constantly to about 18%. How- 
ever, this threshold is important for the tracking detector 
for the following reason: The non-QE processes lead not 
only to single-proton knockout but also to multi-nucleon 
knockout through pion absorption processes and rescat- 
tering. If the proton threshold is zero, these processes 
are not counted because there is more than one proton 
present. Increasing the threshold also increases the prob- 
ability that only one proton is above threshold, in which 
case the event is CCQE-like. Above a certain kinetic en- 
ergy on (« 0.1 GeV), more and more protons are below 
threshold and the ratio decreases again. The dependence 
on the pion kinetic energy threshold is displayed in the 
lower panel. Here the ratio increases because, with in- 
creasing threshold for the outgoing pion, more and more 
non-QE events arc misidentificd as CCQE. 



FIG. 5: (Color onhne) Ratio of the CCl7r+-like to the true 
CCItt"'" cross section as a function of the pion detection kinetic 
energy threshold for CC v^, on ^^C at = 1 GeV. Top 
panel: the solid line is obtained with FSI, the dashed without 
FSI. Lower panel: the different contributions to the full result 
(corresponds to the solid line in the top panel) are shown as 
indicated in the plot. 



We note that a realistic muon-kinetic energy threshold 
of roughly the same magnitude has no visible influence 
on the CCQE-like to true CCQE cross-section ratio since 
the muon kinetic energy is larger in most cases. 



B. CCItt identification 

The CCl7r+ reaction is the second largest cross section 
at the energies of interest in this work, and the major 
background to the CCQE signal channel as we have seen 
in the previous section. 

As in the case of CCQE, the detected CCl7r+ events 
can also be masked by FSI. However, as we will now 
show, the misidentification is minor and independent of 
the detector type — both of our generic detectors identify 
CCl7r+ in the same way. Problematic, however, is the 
low efficiency caused by strong pion-absorption effects. 
The top panel of Fig. [5] shows that already without any 
threshold cuts only 60% of the pions leave the nucleus 
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FIG. 6: (Color online) Distribution of the reconstructed neutrino energy according to Eq. ([Sjl for El'^^^ — 0.5, 0.7, 1.0 and 
1.5 GeV. The reconstructed energy denoted by the dashed lines includes only true CCQE events, while the solid lines are 
obtained by reconstructing the energy with CCQE-like events under Cherenkov assumptions. 



and can be detected.^ Increasing the pion kinetic energy 
threshold decreases clearly the CCl7r+ event rate in the 
detector. In the lower panel, we plot the different contri- 
butions separately and find that the A excitation domi- 
nates. Concluding, Fig.[5]shows that the experiment sees 
only less than 60% of all pions, with that number decreas- 
ing rapidly with increasing pion kinetic energy threshold. 
Therefore, a large part of the total pion yield has to be re- 
constructed. Any data on pion production thus contain 
a major model dependence. This makes it mandatory 
to use state-of-the-art and well-tested descriptions of the 
ttNA dynamics in nuclei. 



V. NEUTRINO ENERGY RECONSTRUCTION 
A. CCQE 

In long-baseline (LBL) experiments, CCQE events are 
commonly used to determine the kinematics. The 
neutrino energy has been reconstructed from QE events 



^ We have normalized the true CCl7r+ to the "no FSI" curve at 
^thros ~ ^^ -^°te also that we use the data of the Argonne bubble 
chamber experiment (ANL) as reference for our elementary pion 
production cross section 



at the MiniBooNE experiment [Hi using 

,ec ^ 2(Ma. - Eb)E^ - {El - 2MnEb + m^) 

2 [{AIn - Eb) - Ef, + \k'\ cos 6^] ' ^ ^ 

with a binding energy correction of Eb = 34MeV and 
the measured muon energy, Ef^, and scattering angle, 6*^. 
The K2K experiment uses the same expression but with 
Eb = [l^. Eq. ([3]) is based on the assumption of 
quasifree kinematics on a nucleon at rest. 

In Fig. [6] we plot the distribution of the reconstructed 
neutrino energy obtained using Eq. ((S]) with Eb = 
34MeV for four fixed El"""^ (0.5, 0.7, 1.0 and 1.5 GeV). 
The dashed lines show the true CCQE events only, the 
solid lines all CCQE-like events (using the Cherenkov 
definition, but without any threshold cuts). Both curves 
show a prominent peak around the real energy which is 
slightly shifted to higher Ef,'"^. This shift is caused by the 
difference between our potential and the specific choice 
of Eb-^ The peak has a width of around 100 MeV full 



^ See Fig. 1111 the dash-dotted and dotted lines there have been 
obtained with Eg = 0. We note that in an event simulation 
the binding energy parameter could be adjusted such that the 
maximum of the reconstructed distribution is at the true energy. 
This is, however, not possible in the actual experiments where 
the true energy is not known. 
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FIG. 7: (Color online) Distribution of the reconstructed neu- 
trino energy according to Eq. ([3]) for _E^"'*' = 1 GeV. Using 
only true CCQE events for the reconstruction leads to the 
dashed line. Including CCQE-like events (Cherenkov defini- 
tion) with various charged pion detection thresholds, one ob- 
tains the solid line (no pion threshold), the dash-dotted line 
(100 MeV), and the dotted line (300 MeV). 



TABLE I: Expected value, E = dEl"" and stan- 

dard deviation, S = (^J^ dEl'"' (g^-g)" ^''^^ foj. ^he 
distributions shown in Fig. [6] and Fig. [S] 



g^"' [GeV] E [GeV] S [GeV] 
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width at half maximum (FWHM). This broadening is 
entirely caused by the Fermi motion of the nucleons — 
Eq. ([3]) assumes nucleons are rest. 

While the distribution of the reconstructed energy for 
the true CCQE events is symmetric around the peak, this 
is not the case for the CCQE-like distribution. The recon- 
struction procedure now includes also non-CCQE events. 
However, Eq. ([3]) is entirely based on the muon kinemat- 
ics and, in the case of A-induced non-CCQE events, more 
transferred energy is needed than for true CCQE, so the 
muon energy is smaller. This lower muon energy leads 
then to the second smaller bump at lower reconstructed 
energies. Thus, the asymmetry is caused by the non- 
CCQE events that are identified as CCQE-like. 

The asymmetry is very sensitive to detection thresh- 
olds, in particular to the kinetic energy threshold for 
charged pions (see Sec. IIIII for the thresholds applied 
in present experiments). We have seen in the previ- 
ous section that increasing this threshold also increases 
the CCQE-like cross section (via the non-CCQE events). 
Thus, a higher threshold leads to a more pronounced sec- 
ond bump, as seen in Fig. [71 

The reconstructed energy under tracking detector as- 
sumptions is plotted in Fig. [H We have seen in the pre- 
vious section that the tracking detector allows the ex- 
traction of a much cleaner CCQE-like sample than the 
Cherenkov detector — almost no fake, i.e., non-CCQE 
events spoil the CCQE-like sample. Consequently, the 
reconstructed distribution is again symmetric, but at the 
cost of a lower detection rate. 

The previous findings, without any threshold cuts. 



arc summarized in Table [If and in Fig. [9l The for- 
mer lists the expected values for the reconstructed en- 
ergy and the standard deviation, while the latter shows 
the probability distribution of the relative discrepancy 
(^roai _ £:roc-)/^rcai f^^. 4 different real energies. We 

note that similar investigations by Blondel et al. 
and Butkevich [isj result in smaller discrepancies. Both 
works consider only CCQE in the initial state, and do 
not include, e.g., A excitation. 

So far, we have discussed the uncertainties in the en- 
ergy reconstruction assuming a fixed, sharp neutrino en- 
ergy. In reality, the energy distribution of the neutrinos 
is broad and thus the question arises how these flux dis- 
tributions are affected by the reconstruction procedure. 
Therefore, we show in Fig. [10] the reconstructed energy 
distribution for the MiniBooNE flux (top panel) and the 
K2K flux (lower panel). Compared to the true CCQE, 
we find an enhancement at low reconstructed energies 
caused by the non-CCQE induced CCQE-like events in 
a Cherenkov-like detector (dashed vs. solid lines, corre- 
sponding to the low-energy bump in Fig.[S]). In a tracking 
detector, the event rates are reduced (dashed vs. dash- 
dotted lines). 



B. CCl7r+ 

The MiniBooNE Collaboration reconstructs the neu- 
trino energy not only using the CCQE sample, but also 



* Note that the standard deviations do not reflect the low-energy 
tails caused by the misidentification of events. 
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FIG. 8: (Color online) Same as Fig. |6l but under tracking detector assumptions for the CCQE-like events (solid lines). 




FIG. 9: (Color online) Normalized distribution of the reconstructed energy vs. the relative discrepancy using Eq. ([3]) for 
El!"^^ = 0.5, 0.7, 1.0, and 1.5 GeV. The dashed lines use true CCQE, the solid ones CCQE-like in a Cherenkov detector and 
the dash-dotted ones CCQE-like events in a tracking detector for the reconstruction. 
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FIG. 10: (Color online) Reconstructed energy distribution for 
the MiniBooNE flux (top panel) and K2K flux (lower panel) 
under different detector assumptions. Eq. Q is used for the 
reconstruction, but with Eb = in the K2K case. 



FIG. 11: (Color online) Distribution of the reconstructed neu- 
trino energy according to Eq. Q for El"^^^ — 1 GeV. Shown 
is the reconstruction based on the CCQE-like sample (before 
and after FSI and Cherenkov assumptions) and based on the 
CCItt"*" sample (before and after FSI). 



f-dEl- and 



TABLE II: Expected value, E 

standard deviation, S = dEl!^'' ^ 

the distributions shown in Fig.[TT]for El"""^ = 1 GeV. 

E [GeV] S [GevP 
from CCl7r+, before FSI 0.94 0.16 (17%) 
from CCl7r+, after FSI 0.95 0.19 (20%) 
CCQE-like, before FSI 0.97 0.13 (14%) 
CCQE-like, after FSI 0.90 0.21 (23%) 



1/2 



for 



using the CCl7r+ sample. Based on the observed muon 
kinematies, treating the interaction as a two-body colli- 
sion and assuming that the target nucleon is at rest inside 
the nucleus, one finds [ll] 

^ _ 1 2MnE^ + Alj - Mj, - ml 
2 -Ef, + cos e^^^El-ml 

where Mjv is the mass of the nucleon, is the mass of 
the muon, 6^, its scattering angle, and its energy. M/ 
is the Breit-Wigner mass of the P33(1232). This formula 
thus assumes that all pious are produced through the 
excitation of the A resonance which is taken to be a state 
of fixed mass, or, in other words, its spectral function is 
taken to be a (5-function. Binding effects are neglected 
here. For Mf = Afjy, this formula agrees with Eq. ([3]) 
for Eb = 0. 

Fig. [TT] shows the reconstructed energy distribution ac- 
cording to Eq. ([3]) for the CCQE-like sample and the 
CCl7r+ sample (before and after FSI). The shape of the 
dash-dotted and the dotted curves have been discussed 
before: Fermi motion broadens the peak and the fake 
CCQE events cause the bump at lower reconstructed en- 



ergies. Also, the reconstructed energy from the pion sam- 
ple is affected by Fermi motion (dashed and solid lines). 
A further broadening comes from the actual shape of the 
A resonance which is taken to be of 5 function-like shape 
in Eq. (|4|). Overall, the reconstructed energy is centered 
around the true energy for both samples, although with a 
slight tendency to lower reconstructed energies. Table HIl 
lists the expected values for the reconstructed energy and 
the standard deviation. Note that the expected value is 
closer to the real energy when using the pion sample and 
that the standard deviation is smaller than in the CCQE- 
like case (calculated here also with Eb = 0). 



VI. q2 reconstruction 

If one assumes a dipole ansatz for the axial form fac- 
tor. Fa, the axial mass. Ma, is the only free parameter 
in the QE nucleon hadronic current (see, e.g., Ref. for 
details; here we use Ma = 1 GeV). Ma affects both the 
absolute value of the cross section and the shape of the 
distribution. Thus, there are two ways of extracting Ma 
experimentally (we assume that the vector form factors 
arc known): (1) Q^-shape-only fit which has the advan- 
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tage that it does not require absolute flux normalization, 
(2) fit to the total cross section. On nuclei, the extraction 
of Ma is much more complicated. Nuclear effects change 
the shape of the distribution and, consequently, the 
extracted AIa depends on the model used to relate mea- 
sured rates on nuclei to nucleonic form factors. Further- 
more, we saw in the previous section that FSI influence 
the CCQE identification. Misidentified events are likely 
to follow a different Q"^ distribution and also affect the 
total cross section as discussed in connection with Fig. [TJ 
Like the neutrino energy, is not an observable — it 
has to be reconstructed from the measured muon prop- 
erties. Using Eq. ([3]), we obtain the reconstructed 
via 



-2£;,(S^-|k'|cos0^) 



(5) 



The neutrino energy itself is reconstructed according to 
Eq. thus Eq. ([5]) is also based on the assumption of 
quasifree kinematics. Fig. [12] shows the CCQE-like 
distribution (solid line) separated into CCQE-induccd 
CCQE-like (dashed line) and fake CCQE (dash-dotted 
line) together with the reconstructed cross section. If 
the background subtraction is perfect (i.e., when the true 
CCQE sample is isolated and only this sample is used to 
reconstruct Q^), then the reconstructed spectrum almost 
reproduces the true spectrum (dashed and double-dashed 
line almost coincide). If background events, namely non- 
QE induced events, are also taken into account for the 
reconstruction ("total reconstructed") then, for the ex- 
treme case that no background at all is subtracted, we 
find an increase at lower Q^, but then it falls off faster 
(dotted vs. solid line). The difference is caused by the 
different muon kinematics of the "fake" events. To con- 
clude, we find that the reconstruction with the simplified 
formulas above turns out to be almost perfect when only 
true CCQE events are taken into account but not if the 
whole CCQE-like sample is used to reconstruct . The 
fake events affect both the height and the slope of the 
distributions and, thus, the extracted Ma values. 



VII. IMPLICATIONS FOR OSCILLATION 
PARAMETER MEASUREMENTS IN v^, 
DISAPPEARANCE EXPERIMENTS 

We close this article with a brief discussion on why 
the exact knowledge of the neutrino energy is of major 
importance. The oscillation probability for the transi- 
tion Va — >■ Vfi is given by (within a simplified two-flavor 
model) 



sin^ 2Q sin^ 



(6) 



where Q is the neutrino mixing angle, Am^ = — m\ is 
the squared mass difference and L the distance between 
source and detector. Consequently, 
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FIG. 12: (Color online) Flux averaged dcr/dQ^ distribu- 
tion of CCQE-like events (solid line) at MiniBooNE condi- 
tions. The dashed line shows the CCQE-like events induced 
by CCQE, the dash-dotted line shows the non-QE induced 
CCQE-like contribution. In addition, the reconstructed spec- 
tra are shown (dotted and double-dashed line). Note that the 
dashed and the double-dashed lines almost overlap. 



The oscillation probability depends directly on the neu- 
trino energy E^^ so measuring the oscillation parameters 
Q and Am^ requires the knowledge of both L and E^. 
In LBL experiments, L and E^, are typically chosen such 
that the detector is placed in the oscillation maximum or 
minimum. Commonly, disappearance experiments mea- 
sure the neutrino flux at the far detector and compare it 
to the one measured at the near detector.^ From the dif- 
ference between both spectra one can determine the oscil- 
lation parameters (compare, e.g., the oscillation analysis 
performed at the K2K experiment [l7|). 

A schematic example is given in Fig.[T3]for 6 « 45° and 
Am^ = 2.5 X 10"'^ eV^ (i.e., the parameters measured 
in disappearance [13 )• It shows the K2K flux with 
L — 250 km. The un-oscillated spectrum is given by the 
dashed line, the survival probability by the dash-dotted 
line, and its convolution by the solid line. An exact re- 
construction of the neutrino energy is thus necessary to 
resolve the oscillated flux, in particular the characteristic 
oscillation dip. 

For further illustration, we show in Fig. [14] the con- 
volution of the oscillation probability with the recon- 
structed K2K energy flux for the three detection scenar- 
ios introduced before: reconstruction using true CCQE, 
CCQE-like (Cherenkov), or CCQE-like (tracking detec- 
tor) events. Clearly visible is the difference between the 
CCQE-like Cherenkov-based reconstruction and the two 
other methods at low neutrino energies around the oscil- 
lation minimum. Note that, in addition, also the recon- 



l-Pn, 



(7) 



This is in contrast to neutrino-appearance experiments which 
measure directly the appearance of a different neutrino flavor in 
the beam. 



11 



K2K flux X (1 -Post;) - 
K2Kflux -- 
(1-Poso) 




FIG. 13: (Color online) Sketch of the influence of neutrino 
oscillations on the neutrino flux (details are given in the text). 
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FIG. 14: (Color online) Reconstructed energy distribution 
for the K2K flux folded with 1 — Pose under different detector 
assumptions (compare to Fig. [10] where the corresponding un- 
oscillated spectra are shown). 



structed energy is uncertain by up to 20% (see discussion 
in Scc.|V|. This effect is important for i/^ disappearance 
searches because it affects the extraction of the oscilla- 
tion parameters obtained from a fit to these distributions 
(see also Ref. [13). 

It is even more relevant when an experiment uses dif- 
ferent detector types for the near and far ones, e.g., the 



K2K experiment uses a tracking detector to measure the 
un-oscillated flux and a Cherenkov detector for the os- 
cillated one. Then one has to extract the oscillation pa- 
rameters from that difference — e.g., by comparing the 
dash-dotted line of Fig. [10|to the sohd hue of Fig. [14| A 
good understanding of the energy reconstruction is thus 
necessary to extract meaningful oscillation results. 



VIII. SUMMARY 

We have applied the GiBUU model to questions rel- 
evant to oscillation experiments. This model provides 
a theory-based and consistent treatment of in-medium 
modifications and final-state interactions. In particular, 
it describes electron- and photon-induced reactions suc- 
cessfully which serve as a benchmark for the neutrino- 
induced processes discussed here. 

The present work addresses the relevance of CC reac- 
tions for neutrino-disappearance experiments. We have 
argued that a correct identification of CCQE events is rel- 
evant for the neutrino energy reconstruction and, thus, 
for the oscillation result. A significant part of CCItt^ 
events is detected as CCQE-like, which is mainly caused 
by the pion absorption in the nucleus. We have found 
that present-day experiments miss the total QE cross sec- 
tion by about 20% and the total pion yield by about 40%. 
These errors have to be corrected for by means of event 
generators so that the final experimental cross sections 
contain a significant model dependence. Furthermore, we 
have investigated the infiuence of these in-medium effects 
on the neutrino energy reconstruction and on the CCQE 
cross section, which is the signal channel in oscillation ex- 
periments stressing the effect of final-state interactions. 

We conclude that any model that aims to describe the 
experimental measurements must — because of the close 
entanglement of CCQE and CCItt"'" on nuclei — describe 
both equally precise, and, in particular, the directly ob- 
servable rates for nucleon knockout and 1tt~^ production. 
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